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Abstract. We evaluate a determinant of generalized Fibonacci numbers, thus pro- 
viding a common generalization of several determinant evaluation results that have 
previously appeared in the literature, all of them extending Cassini's identity for Fi- 
bonacci numbers. 



1. Introduction 

The well-known Fibonacci sequence is given by = fn-i + fn-2 with /o = /i = 1- 
Numerous properties of this sequence are known. We refer the reader to the monograph 
[9] for a wealth of information on this sequence. One of these properties is the so called 
Cassini identity, given by 

fnfn+2 — fn+1 = 

which can be written in matrix form as 

fn fn+ 1 



det 



fn+l fi 



n+2 



-1 



1.1 



Miles [6] introduced k- generalized Fibonacci numbers fn^ by 



(fc) 



EAk) 
J n—i ? 



i=0 



with fi^^ = for every 0<n<k-2, fl^% = 1, and he gave the following generalization 
of (fTTD: 



I Jn Jn+1 ' ' ' Jn+k-1 



det 



n{k) n^K) 

Jn+l J n+2 



(k) 



Ak) 
Jn+k 



\ Ak) Ak) Ak) I 

\Jn+k-l Jn+k ' ' ' J n+2k-2/ 



(2ra+fc)(fc-l) 



1.2) 



More recently, Stakhov [8J has generalized Cassini's identity for sequences of the form 

fn fn—1 fn—p—1- 

Hoggat and Lind [4j consider the so called "dying rabbit problem" , previously intro- 
duced in [T] and studied in [2] or [3], which modifies the original Fibonacci setting by 
letting rabbits die. In previous work by one of the authors [7j, the sequence arising in 
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this setting was studied in detail. For instance, the recurrence relation for this sequence 
depends on two parameters k,£ >2 and is given by 



where Cg , . . . , C^Y£-2 initial values which will be specified below. It was also 
proved that, if ri, . . . ,rfc+£_i are the distinct roots of gk/{x) = x^^^""^ — ^^j-, then the 

k+l-l 

general term of the sequence is given by Cn'^'' = a^rj, with 



i=l 



-1 



^k+e+i-l 




rr - 1 „Yt l>^ r'! - 1 



.«=0 ■ I V I ^1 l=k~l 

Given the previous sequence, for every j > we can define a matrix Aj^k/ by 



A 



■j,k,e 



/ Mk,i) Mk,e) ^(fe/) ^(k,i) \ 

I ^j+e ^j+e+i ■ ■ ■ ^j+k+2e-3 » 

Mk,i) ^(k,i) Mk/) Mk/) 

^j+l '~^j+e+2 ■ ■ ■ ^j+k+2l-2 

. ^{k/) ^{k/) ^(k/) ^(k,e) I 

y^j+k+e-2 ^j+k+2e.-2 ^j+k+2e.-l ■ ■ ■ ^ j+2k+3e.-b- J 

The main goal of this paper will be to find an explicit expression for det(y4j^fc/), thus 
extending (O) and (O). 

2. Extending Cassini's identity 

Before we proceed, we have to fix our initial conditions. In the original setting [7|, 
when we start with a pair of rabbits that become mature I months after their birth 
and die k months after their matureness, the /c + £ — 1 initial conditions are given by 
Ci''') = ■ ■ ■ = Ct? = 1 and Clt''^ = Ct'} + C^n-} fo^ every £<n<k + l-2. Instead, 
in what follows we will consider the following initial conditions: 

= 1, 

^{k/) ^ . . . ^ ^{kl) ^ g 

_ _ p;(fc/) _ , 
•-^fc — ■ ■ ■ — <--fe+<>_2 — 

Note that this change in the initial conditions results only in a shift of indices. Namely, 
if Cn'^'' denotes the original sequence and Cn'^'^ denotes the sequence given by the same 
recurrence relation and these new initial conditions, then for every n > we have 

Mk/) — Piik,e) 

Thus, if Aj^k,e. is the corresponding matrix (defined in the obvious way), we have Aj^k,i = 
Aj+k+i,k/- Hence, we can focus on finding a formula for det{Aj^k/)- 
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First of all, observe that det{Aj^k,i) = det(ylj_i^fc_^) because Aj^^^i is ob- 

tained from Aj_i^k^e by replacing the first row by the sum of the first k rows of the 
matrix, and then permuting the rows so that the first row becomes the last one. If we 
apply this idea repeatedly, we obtain that det{Aj^k^£) = (— (iet(Ao,fc,£)- Hence, 
it is sufficient to compute this latter determinant. 

We shall focus now on computing this determinant, which explicitly is 

/ rik/) ?^(k/) ?^(k/) 7^{k,l) \ 

?;(fc,£) ^{k,e) -^{k.i) ^{k,i) 



det(Ao,fc,£) = det 



\f^{k,i) ^{k,e) ^{k,e) ^{k,£) 

\^k+i-2 ^k+2l-2 ^k+2i-l ■ ■ ■ ^2k+U-b' 



k+l-l 



To do so, recall that we have Cn'^'' = Ost", where the Oj's are given by (11.31) . 



s=l 



We substitute this in the above determinant and use multilinearity in the columns to 
expand it into the sum 



(k+e~i 
n ) ^jfi 



l<i<k+e-l 



^.-1 ^i+e-l i+e . . . i+fc+2£-4 



Now, if in this sum two of the Sj's should equal each other, then the corresponding two 
columns in the determinant would be dependent so that the determinant would vanish. 
We can therefore restrict the sum to permutations of {1, 2, . . . , + £ — 1}. With Sk+e-i 
denoting the set of these permutations, this leads to 



^k+e-i 



<i<fc+£-l V '^(l) ^^(2) 'a(3) ' <T{k+i-l) 



a+e-i \ /k+e-i \ 

i=l / aeSk+e-i V i=2 J - J- 

i=i / <Te5fc+f_i \ j=2 J - j- 

''k+t-1 \ / \ /k+e-1 



n n (-.--^) E (^g--) n -^r" 



j = l / \l<i<j<k+e-l / a€Sk+e-i \ j=2 

I'k+l-l 



i=i / \i<i<j<k+e.-i / 

^k+l-l \ / \ /k+l-l \ 

n% n (-.--^) n-^ 

i=l / \l<i<j<k+£-l / \ 1=1 

X det (^7^-2^+3 ^-k-l+2. ^-k-i+A ... i\ 
l<i<k+t-l ^ * II ) 
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n s n (-.--o(rr^-r7^) n 

j=l J \l<i<j<k+e-l J \ i=l 

xhe^,{r^\...,r^l,_,). (2.1) 

In the last line we have used the following notations and facts: first of all, hm{xi, . . . ,xn) 
denotes the m-th complete homogeneous symmetric function in variables Xi, . . . , x^, 
explicitly given by 

l<ii<---<ini<N 

Furthermore, the Schur function indexed by a partition A = (Ai, . . . , Xn) in the variables 
Xi, . . . ,xn is defined by 



Sx{Xi, . . . ,Xn) 



det x^+^-^j det ix'^^''-' 

^<i,3<N \ J l<i,j<N 



det (x^-') n (^^-^ 

l<i,j<N \ ' J _ 



l<i<j<N 



It is well-known (cf. O p. 41, Eq. (3.4)]) that for A = (m, 0, . . . , 0) the Schur function 
Sx{xi, . . . , xn) reduces to h^ixi, . . . ,xn)- These facts together explain the last line in 
the above computation. 

To proceed further, let us first observe that, by reading off the constant coefficient of 
QkA^)^ we obtain 

k+l-l 

X{r. = {-If"-'. 
1=1 

Furthermore, we have 



^ k+e-l k+£-l 

9kA^) = - = n - = (-1)'+'"' n ^^(1 - ^r'^) 

1=1 i=l 



k+e-l 

i=l 

Hence, we obtain 

oo 



- r. ^x) 



hm{^l ) ■ • ■ ; ^k+e-l)^ 



k+e-l 



m=0 



n (1 - ri^x) 



i=l 

1 



x—1 



1 — X 



= l~x + x'' - x''+^ + ■■■ + 0{x''+^-^). 

In order to evaluate hi^i^r^'^, . . . , r'/^^^^-^), we just have to extract the coefficient of x^~^ 
in the expansion on the right-hand side. This is easy: if ^ — 1 equals a multiple of k 
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then we obtain 1, if £ — 2 equals a multiple of k then we obtain —1, and in all other 
cases we obtain 0. 

We continue evaluating the other factors in fl2.ip . We have 



n i^^ - ri){r-' - rj') = J] 

l<i<j<k+l-l l<i<j<k+£-l 

n 

l<i<j<k+i-l 



n 



i=l 

n I 

l<i<j<k+l-l 



-1) 



k+e 



k+e-i 



Furthermore, we must compute Y\ ^j- To begin with, recall the formula ( 11.3P and 

the fact that C^'''^ = ''^ = ■ ■ ■ = Cf;?.^ = 1 and Cf ''^ = ■ ■ ■ = Cfj? = 0. With this 
in mind, we get 

k+i-l 

n,_-^.k+£+j-l 



i<j<j<fc+^-i 

(-1) 2 TT M \^ r J - 1 



n r+E 



i<j<j<fc+^-i 



Moreover, observe that 



1 >r-^ r", - I 1 r'^ - I I 
h > h 1 = h > + i 

£—2 

= 1 + V ^ + 1 

i=l j 

r\-\ 



Tjir^ - !)■ 

Here, to obtain the second line, we have used the fact that 1 ^ rj is a root of x 

x-l ■ 



k+e-1 



(Tj - U 
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TT f) — ^ 

Now, to conclude we must compute I I — r. To do so, let a; be a primitive 

fj, -1) 

£-tli root of unity. Then 

k+i-l k+l-l I I k+e-1 

n('-,'-i)= nn(o-)=nn 

j=l j=l 1=1 i=l j=l 

^k+e-1 \ n-\ k+t-\ 

iK^^'-i) n n ^'^-^') 

n irj-l)]{-l)^'^'-'^^'-'^[]j9kA^^) 
j=i J \i=i 

Furthermore, gk^u'^) — oj'^'^^'^^~'^) — '^^^ — — Consequently, we have 



Finally observe that 



(i if gcd(£,A;-l) = 1; 



-L-L u'-l 



0, otherwise. 



We can now collect all the work done to obtain the following result. 

Theorem. For all integers k and i with k,i>2, we have 

( tfe-l = akandgcdii,k-l) = l; 
det(Ao,fe/) = < ifi-2 = pk andgcd{i,k-l) = 1; 

[^0, otherwise. 

Corollary. Let kQ,£Q >2 be any integers. Then the following hold: 

i) The sequence {ae}e>2 given by = \ det(Ao,feo,^)| is periodic, and its period is a 
divisor of ko ■ rad(/co — 1). 

ii) The sequence {/3k}k>k given by j3k = \ det(Ao,fc,^o)l eventually zero. 

Proof. i) Clearly gcd(£, /co — 1) > 1 implies that gcd(£ + A;o-rad(A;o — 1), ^o — l) > 1- 
In the same way, if ^ — 1 and £ — 2 are not multiples of fco, then neither are 
£ + fco • rad(/co — 1) — 1 01 1 + kQ- rad(A;o — 1) — 2. Consequently, if q;^ = 0, also 
a^+A;o-rad(fco-i) = as claimed, 
ii) li k > (,0 obviously neither £ — 1 nor £ — 2 can be multiples of k and therefore 
^fe = for every k > io. 

□ 
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